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Abstract 

We introduce a lattice model of dimers with directional interactions as a paradigm of molecular fluids or strongly correlated Cooper 
pairs in electronic systems. The model supports an intermediate phase that is common to both systems. There are two different ideal 
glasses having no moblity since they possess zero entropy. A pairing parameter is introduced to study the geometrical distribution of 
holes in various phases. 



Phillips Q] has provided strong arguments to suggest 
topological similarities between the phase diagrams of 
electronic and molecular systems, and the presence of 
generic intermediate phases with novel properties by con- 
sidering numerous classical molecular systems and quan- 
tum electronic systems. If true, this will be a ma- 
jor progress in our conceptual understanding of glasses 
121 . Unfortunately, these topological similarities cannot be 
demonstrated because of a lack of reliable analytical so- 
lutions. To overcome this, we propose a phenomenologi- 
cal lattice dimer model (on a regular lattice) that describes 
both a molecular and an electronic system. We investigate 
its phase diagram by solving it exactly on a Husimi cactus 
as an approximation of the original lattice. The solution 
on the cactus goes beyond the conventional mean-field ap- 
proximation [3]. Dimers represent the smallest molecules 
that, in the absence of any hole (collectively representing 
hole, solvent, free volume, dopant, impurity or dissociated 
dimers), can get into a unique ordered state not because of 
the regular lattice structure, but because of physical inter- 
actions. On the other hand, monomeric molecules form an 
ordered state merely because of the lattice structure. Thus, 
dimers can be regarded as the simplest model molecules 
whose physics is primarily controlled by interactions and 
not by the underlying lattice, although the latter certainly 
affects it. The extrapolation of the disordered liquid be- 
low the melting temperature Tm gives rise to the super- 
cooled liquid (SCL), which remains disordered at tempera- 
tures where the crystal (CR) is more stable, and represents 
the stationary limit of metastable states (SMS) which can 
be studied by the use of the partition function formalism 
|0| . It is found that the mathematical continuation of the 
SCL entropy S becomes negative (entropy crisis) ol be- 
low a positive Kauzmann temperature Tg. A continuous 
ideal glass transition has to be invoked at Tk to avoid the 
genuine entropy crisis, since states with negative entropy 
cannot be observed in Nature lalZllSll- The SCL freezes 
at Tk and turns into an ideal glass (IG) for T < Tr. The 
glass is inert, independent of the temperature, and has zero 
entropy so that there can be no mobility in it, since mobility 
requires thermodynamically many configurations (positive 
entropy) to change into each other. 

Dimers can also be thought of as representing strongly 



correlated Cooper pairs of electrons in high-T c supercon- 
ductors [9], and have been investigated by various workers 
Since the antiferromagnetic Mott insulator does 
not support the motion of holes created on doping, it is 
thought that there must exist another insulating but disor- 
dered state, the spin-liquid J9J1, which becomes supercon- 
ducting under doping. At higher doping, superconductivity 
disappears, which makes the superconductor an intermedi- 
ate phase 11JI. Similarly, it has become gradually apparent 
in molecular systems (such as Si0 2 , Ge0 2 that are tetrahe- 
drally bonded) that the short-ranged orientational bonding 
plays an important role in giving rise to a liquid-liquid (L- 
L) phase transition 11121 Il3ll to an intermediate disordered 
phase. In ordinary electrolytes, the attractive dipolar in- 
teractions give rise to a similar intermediate phase in the 
form of head-to-tail chains lfl4ll Orientational interactions 
in the dimer model endow the model with enough rich- 
ness to predict a transition that can be identified as a L-L 
transition to an intermediate phase with orientational or- 
der that is different from the crystalline order. The model, 
when applied to Cooper pairs, also gives rise to the colum- 
nar, staggered, and plaquette phases that have been seen 
in the short-r ange d resonance-valance bond model (RVB) 
of Anderson u(\ Hill . The classical dimer model that we 
consider gives rise to both an ordered crystal and two dis- 
ordered ideal glasses at low temperatures; the latter have 
higher energies than the T = crystal. 

Model. We start with a lattice containing N sites, which 
we will take as a square lattice, a bipartite lattice] 7]. A hole 
occupies a lattice site, while a dimer occupies two consec- 
utive sites and the intervening lattice bond. Let E denote 
the energy or the energy eigenvalue, /i the hole chemi- 
cal potential, iVo the number of the holes, and £l(N , E) 
the multiplicity of the state. The rigorous thermodynamic 
treatment of the model is carried out by using the partition 
function Z N given by 

Z N =J2 n (Na,E)r] N ° exp(-/3£), (1) 

where r/ = exp(/?/x) is the hole activity, and (3 = 1/T 
(/cb = 1), T being the temperature. The limit r/ — > rep- 
resents the absence of holes. In the classical treatment, the 
kinetic term is absent; the remaining energy E is taken as a 
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FIG. 1: Columnar (A) and staggered (B) states on Husimi cactus; 
the labels 1,2,.. in (A) and 1-4 in (B), respectively, are introduced 
to capture the ground states; see description in the text. We show 
dimers (thick bonds) in the horizontal directions; they can also 
point in the vertical directions. (C),(D) show the resonance due 
to attractive interaction (dotted line) e p ; (E) shows the axial inter- 
action (dotted line) between colinear dimers. 



sum of various interaction energies. Excess interaction (en- 
ergy e) between a hole and a dimer endpoint is sufficient to 
describe orientation-independent mutual interactions ifl^l . 
which we restrict to nearest-neighbor hole-endpoint pairs 
for simplicity. There is an orientational interaction (en- 
ergy e p ) between a nearest-neighbor pair of two unbonded 
dimer endpoints provided the corresponding dimers are 
parallel (resonance energy, see Figs. 1C,D). If this inter- 
action is attractive (e p < 0), the ground state at T = 
is columnar. If the interaction is repulsive (e p > 0), we 
introduce an additional attractive axial interaction e a (< 0) 
between the endpoints of two colinear dimers (see Fig. IE) 
so that the ground state at T = is staggered. Such 
ground states also play a central role in the short-ranged 
RVB model of high-temperature superconductivity lllOll . 
where the pair of parallel dimers within a square are said to 
resonate. Let N p denote the number of resonating dimer 
pairs, iV a the number of nearest-neighbor axial contacts 
(endpoint-contacts between colinear dimers), and N c the 
number of nearest-neighbor hole-endpoint contacts, so that 
E = 2N p e p + iV a £ a + N c e. We introduce the Boltz- 
mann weights w = exp(—(3e),w p = exp(— (3e p ) and 
w a = exp(— /3e a ) for later convenience. 

The sum in © is over N ,N p , N a and N c , consistent 
with a fixed N. The densities (fik are the limiting ratios 
N k /N as N — > oo; here k = 0,p,a, and c. We always 
take negative /z(< 0) to insure a fully dimer-packed ground 
state at absolute zero. The adimensional free energy uj = 
(1/iV) In Z in the limit N — > oo represents the osmotic 
pressure across a membrane permeable to the dimers, but 
not holes I16I1 . The reduced pressure tt = Pv , where v 



is the volume of a lattice site, is given by ir = Tuj — /i. 
For rj = 0, 7r is not a useful quantity. Therefore, it is 
convenient to introduce the shifted free energy F(T) = 
F(T) — F(0), which remains meaningful for all 7] > 0; 
here F(T) = —Tuj. The equilibrium state must have the 
lowest F(T) among all possible states obtained at given 
T, w p , w a and /i 11611 . The entropy per site s is given by 
s = uj — f3e — <j) In f]- At T = 0, C = since there are 



no holes {q - 
The ratio 
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denotes the number of hole pairs per hole; here, (p Q0 is the 
density of nearest-neighbor contacts between holes. In the 
RVB model, it is a measure of how effective is the pairing 
mechanism for mobile dopants. If all holes were paired as 
isolated pairs, we have D = 1/2. It is commonly believed 
that holes are always paired in the crystal on a square lat- 
tice, but may not be on other lattices lfl7ll . If holes are sepa- 
rately clustered as four sites of a square cell, then D = 1. 
For randomly distributed holes, 4> 00 = 20q on the current 
lattice, and D n = 2(j) . However, the dimer connectivity 
will create correlations even at infinite temperatures, and 
we expect a higher tendency to pair so that D > 2(f> - 
Our results below confirm the additional correlations in the 
disordered phase. It is clear that the values of 4> , and D n 
allow us to draw important conclusions about the way holes 
are distributed in the system. 

Husimi Cactus Solution. In order to solve the model 
exactly, we replace the square lattice by a site-sharing 
Husimi cactus, a recursive lattice, obtained by connect- 
ing two squares at each sites. This is the only approxima- 
tion we make. The cactus can be thought of as a checker- 
board version of the square lattice, representing squares of 
a given color the squares of the other color are miss- 
ing. However, a pair of dimers on the cactus that would 
have belonged to a missing square on the original square 
lattice is counted as a parallel pair on the cactus |4|. The 
model is solved exactly using by now the standard recur- 
sive technique, as discussed elsewhere 10.0, H||. Because 
of its exactness, the calculation (i) respects all local (such 
as gauge) and global symmetries, in contrast to the con- 
ventional mean-field solution which is known to violate lo- 
cal symmetries, and (ii) thermodynamics is never violated 
|2J. The fix-point (FP) solutions of the recursion relations 
represent possible states in the system. The free energy is 
calculated using the method originally proposed by Gujrati 
120- The FP solution that maximizes the osmotic pressure 
represents the three stable states that are found: the ordered 
phase, i.e. the crystal (CR) at low temperatures, the disor- 
dered phase, i.e. the equilibrium liquid (EL) at high tem- 
peratures, and an intermediate phase (IP) with intermediate 
orientational order involved in a liquid-liquid (L-L) transi- 
tion; see below for a complete description. Abandoning 
this maximization principle and continuing the FP solu- 
tions allows us to obtain the stationary metastable states, 
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FIG. 2: (A) e p = -l,e a = 0,(B) e p = l,e a = -1; e = 0.1,77 = 
—0.1 in both cases. 



from which we construct the ideal glassy states 
place the metastable states of negative entropy. 

The exact solution can also be taken as the approximate 
theory for the square lattice. The continuation of EL to 
lower temperatures describes SCL, which as we show here 
exhibits the entropy crisis at T = Tk. The disordered liq- 
uid EL undergoes a transition to IP. The continuation of IP 
to lower temperatures also exhibits an entropy crisis of its 
own at a temperature T = T K that is usually different from 
Tr. We study various contact densities in CR, IP and EL, 
and their continuation. 

Results. We briefly discuss some important results, with 
details found in[4]. For 77 = (fi — > —00), we have ana- 
lytical results, which clearly establish three distinct phases. 
For finite fx, the results are obtained numerically. It should 



be stressed that the existence of IP depends on the choice 
of the energy parameters and the hole chemical potential. 
We note that in Fig. 2, the stable states (EL, IP, and CR) are 
shown by filled symbols (circles, triangles, and diamonds, 
respectively), the realizable metastable states (S > 0) re- 
sulting from EL and IP by the same respective empty sym- 
bols, and unrealizable metastable states (5 < 0) by the 
smaller empty symbols. The point S = denotes the 
Kauzmann point. The extensions corresponding to the un- 
realizable states must be replaced by horizontal lines pass- 
ing through the Kauzmann point (not shown), which rep- 
resent ideal glass states that are frozen in their respective 
metastable state at Kauzmann point. 

Disordered Phase. The disordered phase is charac- 
terized by random orientations of individual dimers, and 
represents EL at high temperatures (T > Tvi) and its ana- 
lytical continuation SCL at low temperatures (Tm > T > 
Tk). The disordered phase free energy, hole density <p , 
and D Q are shown in Figs. 2 A and 2B when the ground 
state is columnar and staggered, respectively. We imme- 
diately note the existence of the entropy crisis in SCL in 
both cases below their respective Tr. The states that rep- 
resent the continuation of SCL over the temperature range 
T < Tk, where it is unphysical due to negative entropy 
must be replaced by an IG at Tc as discussed above. In 
terms of r = [X/w 2 ] 1 ^ 4 , A = 1 + 6u> a u> p + (w a w p ) 2 , 
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2 + w, 



+ u> a + r 2 Wp, 



the analytical expres- 



: are cf)f s = 
,r 2 + 3w a w p + 
We 



sions for thermodynamic quantities for rj = 
(w^r 2 +3w :i w p +w 2 w 2 ) /2r 2 u" , </>p ls = (w ] 
(r 4 + wl)w 2 )/Ar 2 u", andw dls = (1/2) In u" - In 2 
note that the high-temperature limit (w a — ► 1, w p — > 1) of 
Sdis gives (/?Husimi = 1.7071, where ip = exp(2s) is the av- 
erage number of configurations per dimer in the absence of 
interactions. This is closer to the exact value on the square 
lattice, c^exact = 1-7916, [ 18] than that of for Bethe lattice 

VBethe = 1-6875 HI. 

Ordered & Intermediate Phases. The calculation 
for the ordered phase requires separate consideration for 
the two cases. For the attractive case (e p < 0, e a = 0), 
shown in Fig. 2A, the ordered state at T = is the per- 
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1/2. We take dimers in the perfect CR to be 
oriented in the horizontal direction. As T is raised, CR 
begins to distort so that s ord increases, but (/>° rd and <^>° rd 
decrease. As a result, some dimers begin to resonate in 
the vertical directions until resonating dimers in both di- 
rections have equal density at Tp = 1.37, at which point 
CR undergoes a continuous transition into a new state IP 
which is analogous to the plaquette RVB state [20]. This 
phase melts into EL at Tm = 1.58. The continuation of 
EL below Tvi gives rise to SCL, whose entropy vanishes 
at Tk = 0.54. Similarly, the continuation of IP below Tp 
describes its metastable state whose entropy also vanishes 
but at T K = 0.44, thus giving rise to another IG. Both IG's 
have zero entropy, but have different energies and contact 
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densities, so that they represent two different glasses. As 
fx is decreased further from \i = —0.1, the topology of the 
phase diagram does not change; however, the IP entropy 
now has a weaker temperature variation and almost no tem- 
perature dependence in the limit when fj, — > — oo. This 
has the effect that the melting transition at Tvi becomes 
first-order. This limit allows for the analytical treatment, 
which gives <f>? = (3w p + w 2 p )/2X, ^ = 1/4 + 0f/2, 

uj ip = (1/2) In (A) - In y/2, for w a = 1. On the other 
hand, increase of /i(< 0) shifts both Kauzmann tempera- 
tures towards absolute zero, and IP eventually disappears 
0]. Thus, excess holes destroy not only the IG's but also 
IP; the former was also seen in a polymer model recently 
G|. We note that the IG dis and IG IP have almost the same 
hole density as CR. Despite this, D for IG dls has a much 
lower value than for IG IP or CR. It is clear from D = 0.5 
that most holes in the latter are dimerized. In the disor- 
dered phase, even at high temperatures, Dq > 2(f>o , indi- 
cating that there are correlations induced in the hole distri- 
bution due to dimer connectivity. Loosely speaking, holes 
in the disordered phase are mostly isolated, with a tendency 
to dimerize at lower temperatures, which is evident from 
the values of D Q near the Kauzmann temperature where 
4>o = 0; see Fig. 2A. 

The staggered ground state for the repulsive case (e a = 
— l,£ p = 1) is completely frozen for r\ = (fi = — oo), 
and has c/>° rd = 1/2, 0° rd = 0, s ord = 0. It melts at T M 
via a first-order transition into EL. There is a continuous 
L-L transition at Tip in the SCL region where a new IP 
phase emerges. The latter corresponds to configurations 
of dimers preferentially aligned in one direction (vertical 
or horizontal). Its SCL continuation below Tip gives rise 
to the entropy crisis at The IP also exhibits its own 
entropy crisis at The CR state is no longer frozen 
for a finite fj, as shown in Fig. 2B; here, Tm — 1.05, 
Tip = 0.65, T K = 0.41, and T£ =■ 0.56. For a given 
e p , Tip moves towards zero, but Tm decreases to a non- 
zero value Tm = 1.3, as |e a | — > [4]. On the other hand, 
rasing |e a | eventually makes Tp > Tm, hence causing the 
L-L transition to appear in the stable region of the phase di- 
agram [4]. We note that at high temperatures, D > 2(p in 
EL, indicating the presence of correlations in the hole dis- 
tribution due to dimer connectivity. At low temperatures, 
the holes appear dimerized in CR and the two IG's. While 
the SCL hole density is larger then that of IP, < D^ CL 
indicating higher hole pairing tendency in IP. We find that 
Tc as well as Tvi decrease and move to zero as fi (< 0) in- 
creases. Adding too many holes moves melting to T = 0, 
thus destroying the CR state and, therefore, the IG transi- 
tion 0]. 

In conclusion, we have introduced an empirical classi- 
cal model of dimers, and presented its exact calculation 
on a Husimi cactus which serves as an approximate theory 
on a square lattice for molecular and RVB systems. The 
behavior as a function of hole composition can be stud- 



ied by varying the temperature. The calculation provides 
a strong support for the parallel between the two systems 
as observed recently, in particular about the presence of in- 
termediate phases |1J. The empirical model reproduces all 
the known or observed phases in the two systems, which 
justifies its choice. The calculation establishes the exis- 
tence of ideal glasses in the metastable states, and provides 
a qualitative comparison of various stable and metastable 
states. In particular, attention has been paid to study how 
holes are distributed by calculating a pairing parameter D Q 
in various stable and metastable states, which is another 
novel feature of this work. Our work also suggests that the 
ideal glass related to IP may be a strong candidate for the 
spin liquid in the electronic system. We also believe that 
there is a quantum analog of IP that could be relevant in 
high-Tc superconductors. 
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